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A thermodynamic T-matrix approach for elastic 2-body interactions is employed to calculate 
spectral functions of open and hidden heavy-quark systems in the Quark-Gluon Plasma. This 
enables the evaluation of quarkonium bound-state properties and heavy-quark diffusion on a common 
basis and thus to obtain mutual constraints. The two-body interaction kernel is approximated 
within a potential picture for spacelike momentum transfers. An effective field-theoretical model 
combining color-Coulomb and confining terms is implemented with relativistic corrections and for 
different color channels. Four pertinent model parameters, characterizing the coupling strengths 
and screening, are adjusted to reproduce the color-average heavy-quark free energy as computed in 
thermal lattice QCD. The approach is tested against vacuum spectroscopy in the open (D, B) and 
hidden (^I* and T) flavor sectors, as well as in the high-energy limit of elastic perturbative QCD 
scattering. Theoretical uncertainties in the static reduction scheme of the 4-dimensional Bethe- 
Salpeter equation are elucidated. The quarkonium spectral functions are used to calculate Euclidean 
correlators which are discussed in light of lattice QCD results, while heavy-quark relaxation rates 
and diffusion coefficients are extracted utilizing a Fokker-Planck equation. 
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I. INTRODUCTION 

Hadrons containing heavy quarks arc widely used to 
deduce basic properties of the strong force as given by 
Quantum Chromodynaniics (QCD) [l|. Heavy-quark 
(HQ) systems are also valuable for studying hot and 
dense matter, in particular for temperatures and quark 
chemical potentials vifhich are parametrically small com- 
pared to the HQ mass, T,pLq <^ mq (Q = c,b). Under 
these conditions, which are believed to encompass phase 
changes of the medium, HQ momenta {pq ^ mqT) are 
large relative to those of the (light) partons [p^ ~ T^) 
constituting the heat bath. This leads to simplifications 
in the theoretical description which facilitate the task of 
studying the medium. For example, the prevalence of 
elastic interactions with spacelike momentum transfers 
suggest that a potential picture for HQ interactions - 
successful in vacuum - may remain valid in the medium. 

The modifications of heavy quarkonia (charmonium 
and bottomonium) in the medium are believed to re- 
veal quark deconfinement in the Quark-Gluon Plasma 
(QGP) as produced in ultrarelativistic heavy-ion colli- 
sions (URHICs), cf. Refs. @t41 for recent reviews. In 
addition, open heavy-flavor particles are being utilized 
to extract transport properties of the medium formed at 
the Relativistic Heavy-Ion Collider (RHIC), by comput- 
ing their diffusion coefficient in the QGP, see, e.g., Ref. [5[ 
for a recent review. In Ref. Q it was suggested that 
the (in-medium) forces responsible for heavy-quarkonium 
binding may be closely related to those governing the 
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diffusion of an individual heavy quark in the QGP. The 
basis for this idea is that the exchanged 4-momentum, 
k = {kQ,k), for the elastic scattering of a (slow) heavy 
quark is essentially "static", i.e., the energy transfer is 
parametrically smaller than the 3-momentum transfer, 
fco — fc^/277iQ <$i \k I, for both an individual heavy quark 
and a quarkonium bound-state. A possible consequence 
of such a connection could be that a "strong" interaction 
in the medium, which binds charmonium states above the 
critical temperature, at the same time leads to strong cor- 
relations in the heavy-light sector, which accelerate HQ 
thermalization compared to perturbative scattering [6|. 
The purpose of the present paper is to set up and carry 
out a framework which enables a systematic investigation 
of this connection. This requires to evaluate in-medium 
bound and scattering states on equal footing, which will 
be realized by employing a thermodynamic T-matrix ap- 
proach [3, Q- To improve the reliability of this frame- 
work, several constraints will be elaborated: input po- 
tentials will be adjusted to reproduce the HQ free energy 
computed in lattice QCD (IQCD) in vacuum and at fi- 
nite temperature, empirical vacuum spectroscopy in the 
bound-state regime and perturbative QCD in the high- 
energy scattering limit will be checked, and euclidcan 
correlators for heavy quarkonia in medium will be tested 
with IQCD results. 

In the vacuum, the description of heavy quarkonia 
within a potential framework can be made rigorous as 
an effective field theory of QCD with heavy quarks, so- 
called potential non-relativistic QCD (pNRQCD) [l|. In 
a hot and dense medium, however, additional scales enter 
the problem (e.g., temperature, T, and Debye screening 
mass, mi)) rendering the extension of the potential con- 
cept more involved, especially in a strongly interacting 
system where it is difficult to establish scale hierarchies 



(for perturbative treatments based on T ^ mo ~ gT, 
see, e.g., Refs. j9l-[ll|). On the other hand, nonperturba- 
tive information on the HQ interaction over a wide range 
of temperatures is available from thermal IQCD in terms 
of the free energy, F, of a static QQ pair (strictly speak- 
ing, the difference in free energy of the system with and 
without the HQ pair) [l3 - [l4| . In practical approaches, 
the color-singlet free energy, Fi (or the pertinent inter- 
nal energy, Ui — Fi — T dFj&T) has been utilized as 
a potential in Schrodinger |15l - [l8| and T-matrix 0, Q 
equations, and the resulting spectral functions have been 
checked against IQCD results for euclidean correlation 
functions. While these approaches suggest that the po- 
tential model provides a viable tool at finite tempera- 
ture, several open issues remain, e.g.: (i) the use of free 
or internal energy (or even combinations thereof |l9[); 
(ii) the gauge dependence of the color-singlet free en- 
ergy [20[; (iii) microscopic insights into the screening 
mechanisms (e.g., color-Coulomb vs. confining forces); 
(iv) corrections to (or even validity of) the potential 
ansatz. In the present paper we do not offer new prin- 
ciple insights on item (i). To cover the uncertainty as- 
sociated with this problem, our calculations will be car- 
ried out for both free and internal energies which are 
believed to bracket the limiting cases within their in- 
terpretation as a finite-temperature HQ potential. To 
address items (ii) and (iii) we adopt a recently proposed 
field-theoretic ansatz [2l|,l23 to describe HQ free energies 
using a screened Coulomb plus "confining" gluon propa- 
gator. These propagators require four input parameters 
(coupling strengths and screening masses) which are ad- 
justed to reproduce gauge-invariant color-average free en- 
ergies from IQCD. Color projections are extracted within 
the model and utilized to compute color-singlet quarko- 
nium and heavy-light spectral functions, as well as col- 
ored correlations which contribute to the transport of 
heavy quarks in the QGP. Special care is taken of rel- 
ativistic effects - especially for light quarks - which in 
our framework is possible once the vector and/or scalar 
nature of the Coulomb and confining force is specified. 
For example, the relativistic Breit correction known from 
electrodynamics [2j] naturally emerges as a relativistic 
effect in the Coulomb potential. We will furthermore 
check the static approximation underlying the poten- 
tial picture, by comparing different versions of the 3- 
dimensional reduction scheme to obtain the T-matrix 
from an underlying Bethe-Salpeter equation. 



This article is organized as follows. In Sec. [IT] we set 
up the microscopic model used to fit IQCD free energies. 
In Sec. mil we recollect main elements of the thermody- 
namic T-matrix formalism, including relativistic correc- 
tions. Section HV] is devoted to the discussion of our nu- 
merical results for quarkonium and heavy-light spectral 
functions and their applications to euclidean correlators 
and HQ relaxation times, respectively. We conclude in 
Sec.lVl 



II. MICROSCOPIC MODEL FOR THE 
HEAVY-QUARK POTENTIAL 

The recent revival of potential models to describe 
heavy quarkonia in medium has been largely driven by 
the prospect of a paremeter-free input via static HQ free 
energies computed in thermal IQCD. However, functional 
fits to the lattice "data" usually do not offer much insight 
into the physical mechanisms underlying the medium ef- 
fects in the potential, nor do they allow to define ver- 
tex structures of the interaction which become important 
at higher energies and/or in different color channels. In 
addition, it is desirable to base the starting point on a 
gauge- invariant quantity, i.e., the color- average free en- 
ergy^. To this end, we adopt the microscopic model de- 
veloped by Megias et al. [2l|, [22| which we briefiy review 
in the following and then fit to lattice data. The key 
idea underlying this model is that the HQ free energy 
can be accounted for by a nonperturbative ansatz for 
the gluon propagator giving rise to a string-like confin- 
ing term in coordinate space, plus a "standard" pertur- 
bative term corresponding to a screened color-Coulomb 
interaction |2l| . 
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to be understood in static gauge and within dimensional 
reduction (tod, rho'- screening masses). The leading 
nonperturbative effect in the gluon propagator is asso- 
ciated with a dimension-2 gluon condensate dictated by 
dimensional considerations. 
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with ttiq a "glueball" mass. A priori, a dimension-2 
gluon condensate is gauge-dependent and as such a some- 
what controversial quantity. Since the gluon propagator 
is a gauge-dependent quantity the appearance of gauge- 
dependent terms is inevitable. However, it has been ar- 
gued by several authors J24l - l28t that a dimension-2 con- 
densate encodes nontrivial gauge-invariant information, 
e.g., topological configurations associated with magnetic 
monopoles giving rise to a static confining force (which is 
precisely the effect modeled in the present context). Evi- 
dence for a dimension-2 condensate contribution has also 
been found in QCD sumrulcs (see Rcf. [29| and references 
therein) . 

To establish the connection to the HQ free energy 
(given by a correlator of Polyakov loops, fi), one starts 



^ Since the model adopted here and the lattice calculations use 
different gauges (static vs. Coulomb), this is the only meaningful 
way to extract parameters. 



from its perturbative expression at finite temperature in 
the color-singlet channel j2l| . 
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with the perturbative gluon propagator 
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The separation-independent term, {A^^{x)), in Eq. ^ 
plays the role of a sclfenergy of an individual heavy quark. 
The main assumption consists now of augmenting the 
perturbative propagator by the nonperturbative one as 
given by Eq. ([T]). Assuming further that the same func- 
tional dependence as in Eq. ([3]) holds in the nonperturba- 
tive case we are lead to the following form of the singlet 
free energy {Nc = 3): 
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In Ref. [2l[ this approach has been applied to study the 
Wilson loop and HQ free energy in quenched QCD at 
finite temperature, and efficiently reproduces pertinent 
IQCD data. As an extension of this treatment we al- 
low for different screening masses in the perturbative and 
nonperturbative parts of the gluon propagator which im- 
proves the precision in our fits to unquenched IQCD data. 
As indicated above we fit the color-average free energy. 
Since we aim at a parametrization over a large ran ge i n 
distance and temperature we employ the definition |3C| 
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without further approximations which automatically en- 
sures the correct behavior in the limits rT >> 1, 
where the potential is dominated by 2-gluon exchange, 
Fav.{r,T) « (i^i(r,r))2, and rT « 1, where we have 
Fav{r,T) « Fi(V,T) + Tln(9) [U. The evaluation of 
Eq. ([6]) requires the color-octet and -singlet free ener- 
gies (in addition we will use the extracted potentials in 
the sextet and anti-triplet channels for the calculation 
of HQ relaxation times). In previous works [g, Q these 
potentials have been approximated by Casimir scaling of 
the singlet potential. While this is a good approxima- 
tion for the short-r ang e (perturbative/Coulombic) part 
of the potential [3ll |32| | , it presumably does not apply to 
the confining part which rather appears to be universal, 
including its long-distance limit [30, 33, i3J]^. We there- 
fore apply Casimir scaling only to the Coulomb part of 



the model gluon propagator while the string part is as- 
sumed to be color-blind, i.e., the same in all color chan- 
nels^. This is also compatible with the interpretation 
of the long-distance limit as an individual HQ selfenergy, 
as discussed below. The coordinate-space potential in a 
color-channel a takes the form 
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with the Coulomb part F^ , the nonperturbative string 
part F"^, and a r-independent part, Foo{T), which will be 
associated with a (real part of the) HQ "selfenergy". Eg, 
below. Similar analytic forms of the potential have been 
used for fits to the color-singlet potential in Ref. jlSj . 
Let us examine two limits of this expression. First, for 
T — > both screening masses should vanish while the 
condensate characterized by rac remains finite [2l|; one 
obtains 1211 
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which recovers the Cornell potential in the color-singlet 
channel and yields a universal string tension in all color 
channels, consistent with Refs. [30, ISJ, [SJ] . Second, for 
r — >■ oo at finite T one has 
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which is independent of a, consistent with IQCD data 
in Refs. [IJj IS^- When additionally taking the zero- 
temperature limit of Foo, it diverges since ma > and 
rriD — >■ 0. This is, of course, expected in quenched QCD 
but needs to be amended in the presence of light quarks, 
to simulate "string breaking" . We implement this by en- 
forcing a flat potential above a string-breaking scale of 
about r ~ 1.1 fm where the potential has reached about 
1.1 GeV. We now also see that the cancellation between 
the leading piece of the second term and the constant 
third term in the parentheses of Eq. (O in the rhu — > 
limit only works for all color channels if the string term 
is color-blind (i.e., has no Casimir scaling). If, on the 
other hand, both terms proportional to rriQ are subject 
to Casimir scaling, it would imply that the r ^ oo limit 



^ Calculations reported in Ref. [3511 come to a different conclu- 



sion, possibly because they cover smaller distances than those in 
Refs. |30l . 13411 . In Ref. |30|| potential problems with the compu- 
tation of the octet free energy on the lattice below Tc hayc been 
pointed out. 
^ Such a decomposition is not possible in functional parametriza- 
tions of the potential. 
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FIG. 1: (Color online) Color-average HQ free energies obtained in thermal IQCD computations by Kaczmarek et al. (left 
panel ) [12, l3g| and Petreczky et al. (right panel) [13 . ISTI ]. compared to our fits using the microscopic model suggested in 
Refs. Ulnif?^ 
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FIG. 2: (Color online) Temperature dependence of the parameters deduced from our fit to the color-average free energies from 
the two IQCD calculations displayed in Fig. [l] Left panel: strong coupling constant as; Right panel: screening masses of the 
color-Coulomb and confining parts, and dimension-2 condensate "glueball mass", mo- 



(i.e., the single HQ selfenergy) picks up a strong depen- 
dence on the color orientation of the quark, which is not 
natural. 

Once the temperature dependence of the parameters 
mc, tud and mjj, as well as of the strong coupling, as, is 
determined through fits of the color-average free energy, 
corresponding expressions for the internal energy, U, can 
be computed via 



U{r,T) = F{r,T)~T—F{r,T) 



(10) 



and projected into different color-channels, a. 

It is currently an unsettled question whether the free 
or internal energy (or a linear combination thereof) is a 
more suitable quantity to be utilized as a static two-body 
potential in a Schrodinger and/or scattering equation. 
We recall that the quantity, F{r,T), computed in ther- 
mal IQCD, is the difference between the free energies of 
the thermal system containing a static QQ pair and the 



thermal system without the pair. In Ref. [38| it has been 
argued that the pertinent difference in internal energies, 
U{r, T), recovers the thermal expectation value of the po- 
tential energy between the static Q and Q charges. This 
suggests U{r) as the appropriate in- medium two-body 
potential. Such a potential would by construction be a 
real quantity and thus a natural starting point to be uni- 
tarized in a scattering equation, generating appropriate 
on-shell cuts in the intermediate state through imaginary 
parts in the scattering amplitude. In Ref. [39[ it has been 
argued that the relevance of F vs. U depends on the in- 
terplay of the thermal relaxation time in the heat bath 
and the interaction time of the Q and Q. If the for- 
mer is much shorter than the latter, the QQ motion will 
be adiabatic and the free energy should be used; on the 
other hand, for very short interaction times (e.g., for a 
broad resonance or high-energy scattering), the internal 
energy should be more suitable. Another approach to the 
problem has been pursued by using effective field theo- 
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FIG. 3: (Color online) Temperature dependence of the diinension-2 condensate C2, Eq. (|14[l. (left panel) and of the infinite- 
distance limit of the free and internal energies for the two IQCD computations displayed in Fig. [l] 





1 

0.5 



-0.5 

-1 


!■ 


'.'II ;^H> 






\..^:^: : 




> 

CD 


'-- 






S 


'-.^ ^\ -—. 


LL 


/ 8 

/ 6 

/ 3 






il 




3 


0.2 0.4 0.6 0.8 1 
r[fm] 


1.2 



> 

CD 
O 




0.2 0.4 0.6 0.8 

r[fm] 

FIG. 4: (Color online) Vacuum HQ potentials in different color configurations of QQ (singlet and octet) and QQ (triplet and 
sextet) channels, extracted from model-fits to color-average QQ free energies in IQCD for Nf = 2-1-1 \VA. ISg. l4l| ("potential 
1", left panel) and for iV/ = 3 [H, H^] ("potential 2", right panel). 



ries at finite temperature [9|, lUj , by combining HQ and 
perturbative hierarchies (HQ speed v <^\ and gT <^T). 
These studies recover the color-Coulomb part of the in- 
teraction and suggest that the free energy plays the role 
of a potential in a Schrodinger equation. In addition, an 
imaginary part of the effective potential has been iden- 
tified [lO|, |40|- We emphasize, however, that within a 
T-matrix approach, imaginary parts are generated via 
the unitarization procedure in the intermediate propa- 
gator and thus imaginary parts in the potential should 
be implemented via suitable cuts in a coupled-channel 
treatment. To account for the present uncertainty in the 
identification of an irreducible 2-body HQ potential we 
will show numerical results for both F and U as driv- 
ing kernel in the T-matrix equation. More precisely, we 
utilize the subtracted quantity 

Vair;T) ^ Xa{r,T)-X{r ^ 00, T) , X ^ F or U (11) 

to ensure convergence of the Fourier transform. These 
choices are believed to bracket the range of interaction 



strength in the HQ sector. 

Besides the potential, another important ingredient to 
the two-body scattering equation arc the in-medium self- 
energies of the heavy quark and antiquark, which we 
treat symmetrically as appropriate for a hot medium 
at vanishing baryon chemical potential. Starting from 
a bare quark of mass itiq, we interpret the potential 
value at infinite distance, X{r ^ oo,T) in Eq. pTj) . as a 
temperature-dependent "mean-field" contribution to the 
HQ masses, i.e., as a real part of the selfenergy, 

mQ=m°Q + i:^{T) , I]|(r) = X(r^cx),T)/2. (12) 

This interpretation follows from the picture that at infi- 
nite distance the quarks have become independent from 
each other which is supported by IQCD results as dis- 
cussed above. In addition, we will investigate the effects 
of imaginary parts of the HQ selfenergy, 

rQ = -2ImI]Q, (13) 

associated with HQ rescattering in the heat bath. This 
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quantity has been estimated to be about ■ 
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Let us now turn to the fit of the potentials to recent 
IQCD data for the color-average free energy. To obtain 
an indication of the systematic uncertainty underlying 
different IQCD inputs we will carry out all our calcula- 
tions for iV/=2-f 1-fiavor QCD [ll,|3i,|4l| and for A^/=3- 



flavor QCD [13|, l37| (the latter input has been used in the 
color-singlet channel in a previous T-matrix study [7l]), 
which wc refer to in the following as potential 1 and 
2, respectively. The underlying (pscudo-) critical tem- 
peratures in these IQCD calculations have been quoted 
as rc=196MeV (potential 1) and 190 MeV (potential 2). 
Fig. [T] shows the pertinent color-average free energies to- 



gcthcr with our fits which have been performed down to 
temperatures of ca. 0.8 Tj, (not aU are shown in tlic plot). 
The agreement with each data set is fair and supports 
the adequacy of the underlying model. The tempera- 
ture dependence of the four fit parameters is displayed in 
Fig. [2] The variation of all parameters between the two 
potentials is rather small. The strong coupling constant, 
as{T), depends weakly on T. To suppress fluctuations 
in an unconstrained fit, we have for simplicity adopted a 
linear ansatz^. The screening masses exhibit an appre- 
ciable increase with temperature reminiscent of the linear 
T-dependence one expects from leading-order perturba- 
tion theory, m^'"'* = {1 + Nf/Gy/^gT. Compared to the 
perturbative result, the coefficients in our fits are smaller 
for the Coulomb part (ttid) and larger for the "string" 
part {rhu), suggesting a stronger screening of the confin- 
ing part (which primarily acts at larger distances). We 
have tried fits enforcing the condition mu = rhu but 
could not obtain satisfactory accuracy without introduc- 
ing unnaturally large variations of the parameters. This 
might support the assertion of differences in the screening 
process for the perturbative and nonperturbative compo- 
nents of the free energy. On the other hand, the variation 
of the glueball mass with T is weak. Recalling the explicit 
relation to the dimension-2 condensate UM 
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we find that a constant of about C2 — 0.8 GcV^ above Tc 
(see left panel of Fig. [3]) is well compatible with our fit, 
as also found previously in Ref. l2l| and in analyses of 
the gluon propagator, three-gluon vertex or quark propa- 
gator (see Ref. [23 and references therein). On the other 
hand, we have verified that the 20% decrease of C2(T) 
across Tc is a robust feature within "reasonable" varia- 
tions of the other fit parameters; e.g., when imposing an 
overall T- independent value for C2, we could not repro- 
duce the IQCD values of F^ over the entire temperature 
range without "unnaturally" large variations in the other 
fit parameters. It is tempting to speculate that the 20% 
drop in C2 (T) across Tc is related to a similar drop found 
in the magnetic-monopole density in SU(2) gluodynam- 
ics in Ref. |25|, where qualitative arguments have been 
put forward that an Af^ condensate could be connected 
to the monopole density. The temperature window over 
which the variation of C2{T) occurs basically conicides 
with where rapid changes in the infinite-distance values 
Fao and Uoo are observed, cf. right panel of Fig. [3] In 



** Without this constraint, the fits tend to generate what we beUeve 
are artificially large variations in the parameters, mainly caused 
by varying ranges in r as covered by the IQCD data at different 
temperatures. This is particularly evident when fitting to the 
color-singlet free energy and removing some of the small-distance 
points. As a general guiding principle in the fits we tried to utilize 
redundancies in parameter choices to obtain smooth variations 
with T. 



the zero-temperature limit, assuming that both screening 
masses go to zero, a value of m^ ~ 1 GeV is needed to 
reproduce the vacuum string tension of ^/ct = 0.465 GeV 
found in IQCD |36l . |37| . in connection with strong cou- 
plings of as=0.285 and 0.33 for potentials 1 and 2, re- 
spectively. All of these values are close to the fitted ones 
at the lowest temperature. 

In Fig. |4]we summarize the results for the vacuum po- 
tentials in the four different color channels which can be 
formed in 2-body QQ and QQ systems (recall that in 
vacuum the entropy-term vanishes and thus Fa = Ua] 
also Fav = Fi from Eq. ^). The color blindness of the 
string term produces a long-range attraction in all chan- 
nels (which will support colored bound states in vacuum 
as discussed in Sec. lIV A)) . The potentials emerging from 
the model-fit at finite T are collected in Fig. \5\ One 
clearly recognizes the "melting" of the string term with 
increasing T. The singlet (meson) and antitriplet (di- 
quark) potentials remain attractive at all distances. For 
the octet and sextet potential some residual attraction 
from the string term persist at lower temperatures (espe- 
cially in U), preserving a shallow dip structure for quark 
separations of around 0.1-0.2 fm. This behavior has also 
be seen on the lattice [Tj] and is obviously incompatible 
with a Casimir scaling of the string term. 



III. THERMODYNAMIC T-MATRIX AND 
OBSERVABLES 



A. Reduction Scheme and Relativistic Corrections 



The above constructed in-medium potentials are now 
implemented into a thermodynamic T-matrix approach. 
The latter follows from a 3-dimcnsional (3D) reduction 
of the Bethe-Salpeter Equation in ladder approxima- 
tion [4a - [44| . Heavy-quark systems are particularly suit- 
able for this reduction as their energy transfer is para- 
metrically suppressed compared to the momentum trans- 
fer. Even for heavy-light systems the on-shell condition 
on the heavy quark suppresses the energy transfer rela- 
tive to the 3-momentum transfer. Note that a 4D treat- 
ment cannot improve the accuracy as long as the input is 
based on static potentials. However, relativistic correc- 
tions, as well as different reduction schemes, should and 
will be addressed below. The former are necessary to 
ensure consistency for relativistic energies (and, in fact, 
establish "minimal" Poincare invariance of the potential 
approach [4g (see also [4y]))' while the latter give an 
indication of uncertainties inherent in the static approx- 
imation. The 3D integral equation for the T-matrix can 
be further simplified by applying a partial-wave decom- 



position which leads to the following ID equation, 



Ti^aiE; q', q) = Vi,a{q\ q) + - dk P Vi^q', k) 



xGi2(-B; k) Ti,a{E; k, q) [1 - ni.(wi(fc)) - nF(c^2(fc))] , 

(15) 

for the T-matrix Ti^a in a given color channel (a) and par- 
tial wave (l); np is the Fermi-Dirac distribution, q = \q\, 
q' = \q' \ and fc = |A: | denote the relative 3-monientuni 
moduli of the initial, final and intermediate 2-particle 
state, respectively, and u;i{k) = {mf + fc^)^/^ the single- 
quark energies. Equation (|15p encompasses both the 
heavy-light {1 ~ Q, 2 = q) and quarkonium (1,2 — Q) 
channels where either particle can be a quark or an anti- 
quark. The precise form of the two-particle propagator, 
Gi2, depends on the reduction scheme [43, |43: fo^' which 
we will investigate two well-established options, namely 
the Blancenbecler-Sugar (BbS) [42| and the Thompson 
(Th) m scheme. 



GjnE;q) 

Gf^HE;q) 
m{q) 



m{q) 



E-uji{q)-ujQ{q)-Y.q~j:^ 



Q 



2m{q){uj^{q)+ujQ{ci)) 
E^-{io^{q)+uQ{q) + ^, + Y.Qf 

ioi{q)ojQ{q) ■ 



(16) 



The main difference between both schemes is that the de- 
pendence on the total energy, E, is quadratic for the BbS 
propagator but linear for the Thompson version. The 
form of the propagators in Eqs. ([1^ further implies that 
both quarks remain good quasiparticles in the medium, 
i.e., their widths Fg ^ are small compared to their mass. 
We use a minimal width of Fg.g = 20MeV to facilitate 
numerical stability, unless otherwise stated. The incor- 
poration of microscopic quark spectral functions will be 
carried out in an upcoming study. Once the potential is 
specified, the T-matrix equation (1151) is solved using the 
algorithm of Haftel and Tabakin [431 ^^ i^- previous works 
in our context [3, Q . 

It remains to specify how we implement the coordinate- 
space potential as extracted from IQCD in the previous 
section. We start by performing the Fourier transform 
and partial- wave expansion according to 



V^'fiq^q) = 



Xq'q 



'^'"^^'''''Piixq,q)V,^/'ir)e^'^'^-'^y 



Stt 



q -q 

^/WW\ 



(17) 



with the usual Legendre Polynomials Pi . Since the string 
part of the potential {V^) is primarily active at long dis- 
tances, i.e., at low momenta and thus in the nonrela- 
tivistic regime, no further amendments are applied to it. 
However, for the Coulomb part {V'~^), which dominates 
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FIG. 6: One-gluon exchange diagram for quark-quark scatter- 
ing; in the center of mass system, the relative 4-momentum 
in the incoming (outgoing) state is q = (pi — P2)/2 {q' = 
(p1-P2)/2). 



at small distances (and thus at relatively large momen- 
tum transfers), several corrections are in order. To infer 
relativistic effects, let us back up to the definition of the 
relativistic Coulomb potential given by the perturbative 
one-gluon exchange diagram in Fig. |6l Suppressing the 
color structure, the Born amplitude (potential) is given 
by 



V 



Anas 



u{p[)j''u{pi)u{p'2)jf,u{p2) (18) 



•■D 



Yukawa Spinor-contraction 

where we have included a Debye mass as the leading tem- 
perature correction in the gluon propagator. In addition 
to the standard Yukawa propagator (which in the static 
approximation amounts to setting t — > — (<f — q')"^ in 
the center of mass system) we have a contribution from 
the contraction of the spinors with the vertex. At the 
level of the cross section (or amplitude squared) this part 
gives rise to the following factor (with the normalization 
uu = 1), 



\u{p[) 7^ u{pi) u{p'2) 7^ u{p2)\" 
8(2(s- 



,)2+2.sf + f2) 



16 rn^ rriQ 



(19) 



For large s ~ E'^ the terms proportional t are sublead- 
ing and can be dropped so that we can reformulate this 
expression as 



l"(Pi) 7^ u{Pi) u{p'2) 7a* "(P2)|^ 



8(2(s- 



mlf) 



16mi rriQ 



(20) 



= 4 



2 2 



(1 



'Jq^Q 



-r 



The factor in parenthesis is precisely the well-known 
Breit interaction in electrodynamics, corresponding to 
a magnetic current-current interaction of 2 moving 
charges [2j, [S^, [50] , while the first factor "corrects" for 
relativistic kinematics, (the extra factor 4 arises from the 





0.8 




0.7 




0.6 


JH 

^ 


0.5 
0.4 


D 


0.3 




0.2 




0.1 








pQCD 

non rel. 

with corr. 




=1. 



2.5 



3.5 
s^^^ [GeV] 



4.5 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




pQCD 

non rel. 

with corr. 



heavy-heavy 




2.5 



3.5 
s^^^ [GeV] 



4.5 



FIG. 7: (Color online) Left panel: perturbative QCD cross section for quark-quark scattering via i-channel gluon exchange 
using a vector interaction. The full pQCD result (solid line) is compared to the T-matrix result in Born approximation using 
(i) a nonrelativistic treatment of the spinor structure (doted line) and (ii) the correction factors R and B. Right panel: the 
same but for a scalar interaction. The heavy- and light-quark masses are set to 1.7 GeV and 0.4 GeV, respectively, and the 
strong coupling and Debye mass have been fixed at 0.3 and 0.67 GeV, respectively. 



summation over spins and has to be taken out in a spin- 
independent definition of the potential). We therefore 
identify the following factors with which the nonrelativis- 
tic (off-shell) potential, V{q,q'), should be augmented: 



R{q',q) = m(g)-l/2 m(<z')-'/' 
B{q\q) = fe(g)l/2 fe(q')i/2 

b{q) = (1 + 



^q{q)^Q{(i)' 



(21) 
(22) 

(23) 



Note that i?, i? — > 1 for q,q' <C mq^q. For the string 
term, for which we assume a scalar interaction, the spinor 
contraction leads to 



W)u{q)u{p')u{p)\' 

4(4m^ - t){Am^Q 



t) 



16 rrii rriQ 



(24) 



Assuming again that we can drop the terms proportional 
t (relative to ruQ^q), no relativistic correction factor is 
mandated (this refines the procedure adopted in earlier 
works 0,i). 

To check the impact of our corrections we compute 
the cross sections for one-gluon exchange (Fig. ^ for 
quark-quark scattering using the Coulomb term in Born 
approximation including our correction factors. 



da 
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mlm^Q 
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TLa - VP, 



(25) 



La 



and compare it to the exact perturbative QCD (pQCD) 
results in the left panel of Fig. [71 One finds that the 
relativistic correction factors B and R are essential to es- 
tablish consistency with pQCD at high energies; even at 
low energies, the agreement is no worse than 20%, which 
supports the approximation of neglecting t against s and 



niQ^q in the numerator of Eq. (J19p . The factors R and 
B provide a substantial improvement over not includ- 
ing them. Without the R factor, one obtains vanishing 
cross sections at high energy (only half of the correct 
magnitude without the Breit correction); even close to 
threshold the discrepancy to pQCD is larger than with 
the corrections. Also note that without the relativistic 
factors the cross section goes to zero for niq — > 0, which 
becomes an uncontrolled feature in applications to heavy- 
light scattering. In the right panel of Fig. [7] we compare a 
"pQCD" calculation assuming a scalar vertex structure 
to the T-matrix Born results with and without correc- 
tion factors B and R. In this case, it is obviously man- 
dated not to include these factors. As to be expected, 
the nonrelativistic approximation leads to the same re- 
sult irrespective of whether one uses a vector or scalar 
interaction. 

Finally, we account for effects of the running coupling 
constant in the off-shell extrapolation of the potential. 
For on-shell kinematics, q = q' , such effects are effectively 
taken care of by our parametrization of the potential. For 
off-shell scattering, we simulate the running coupling by 
introducing a factor Fi^niq t^ q') < ^ as 



Fruniq',q) = In 



A2 

A2 



/In 



id 



.M2 



A 



21 



A2 



(26) 



with A = 1 GeV and A = 0.2 GeV. 

Putting all corrections together, the final form of the 
potential figuring in the T-matrix equation p5p is 



Vi,a(g',g) 



R{q',q)B{q\q)Frun{q',q)Vi'^,{q',q) 
+Vii{q\q). (27) 



In vacuum the unscreened Coulomb potential exhibits a 
well-known infrared singularity. We tame this by intro- 
ducing a small low-momentum cutoff; we have checked 
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that varying this cutoff has a very small effect on the 
vacuum spectral functions. 



B. Quarkonium Correlators and HQ Diffusion 

The key quantity for computing observables is the on- 
shell T-Matrix, Ti,a{E]q,q), where E = a;i(g) + uj2(q) 
for both in- and outgoing channels. Following Ref. [3], 
the continuation below the 2-particle threshold, i^thr = 
mi + 7712, is carried out for vanishing 3- momentum, 
Ti^a{E;0,0). The mesonic spectral function in a given 
quantum-number channel a is defined as 



aa{E) = --Img^{E), 

TT 



(28) 



where G denotes the correlation function, for which we 
will focus on the case of a heavy quark and antiquark 
(QQ) in the color-singlet channel (a=l) in a QGP of 

It can be written as 



vanishing net baryon charge (fiq) 

g^{E) ^ g"^{E) + AGaiE) 



(29) 



where 

GliE) - iNfN, 



cPk 



GuiE;k) [l-2nF(wQ(fc)] 
(30) 



(27r)3 
xTr TaA+ik)TcA-i-k) 



denotes the noninteracting contribution with the parti- 
cle/antiparticle projectors 



A±(fc; 



2mQ 



(31) 



The Dirac matrices Ta e {IiTmi 75' 7/^75} characterize 
the scalar, vector, pseudoscalar and pseudovector chan- 
nels, respectively (corresponding to Xco, J/'4', Vc and Xci 
in the charmonium sector). In the following we will ne- 
glect effects due to spin-orbit and spin-spin (hyperfine) 
interactions which is justified in the HQ limit. It im- 
plies degeneracy of the S'-wave (^==0) states J/V' and rjc, 
as well as of the P-wave (^=1) Xc states (in the vacuum 
spectrum, this is realized within ^0.1 GeV). Thus, we will 
compute only the vector (Fq, = 7^) and scalar (Fq, = 1) 
channels. The interacting contribution to the correlator 
in Eq. (|29p is given in terms of the off-shell T-matrix as 



AGaiE) = 



NfN, 



dk e Gi2(E; fc) [1 - 2 npiiUQik))] 



X / dk' k' Gi2{E;k')[l - 2nF{ujQ{k'))] 

X [ao{k, fc') To^a{E; fc, fc') + ai(fc, fc') Ti,a{E; fc, fc-')] . 

(32) 

The coefficients ai result from the traces over the spinor 
structure. In line with the above approximation of ne- 
glecting spin-induced interactions, we use a I/itiq ex- 
pansion for these coefficients (which also leads to the de- 
generacy of pseudoscalar- vector and scalar-axialvector). 



One has 



a-o 



ai =0 



for S waves and 

Go — fli 



== -2- 



fcfc' 



(33) 



(34) 



'Q 



for P waves. 

To test or quarkonium spectral functions against IQCD 
correlators, computed with good accuracy in euclidean 
space-time, we need to calculate the euclidean-time cor- 
relator defined as 



G„(T,r) = JdEa^{E,T)ICiT,E,T) , 



(35) 



/C(T,i?,T) 



cosh [E (r - /3/2)] 
sinh [EI3/2] 



The use of a constant width in the two-particle prop- 
agator, Eqs. (|16p . implies an non- vanishing value for 
aa{E — > 0). This induces an artificial singularity in 
the euclidean correlator since the temperature kernel, 
K., diverges in the zero-energy limit. This is an arti- 
fact of the quasiparticle approximation that can in prin- 
ciple be cured by employing a microscopic calculation 
of the single-particle selfenergies leading to the proper 
limit, cra{E — s- 0) ^- -E, for the retarded meson spec- 
tral function. We defer this study to future work and 
evade the singularity problem by imposing a cutoff be- 
low which we set the imaginary part of the propagator 
to zero, i^cut = 2(8) GeV for charmonia (bottomonia); 
there is very little sensitivity to our correlator results 
when decreasing E'cut by up to a factor of 2. 

An additional subtlety in the comparison of model 
spectral functions to IQCD euclidean correlators [5l|, [54I 
is the presence of so-called zero-mode (ZM) contribu- 
tions. These may be pictured as changing the time di- 
rection of a HQ line and thus represent HQ scattering 
processes including Landau damping (rather than QQ 
propagation). It turns out that the pseudoscalar quarko- 
nium channel does not pick up the ZM contribution. To 
avoid extra uncertainties due to the latter, we therefore 
restrict our comparisons to euclidean IQCD correlators 
to this channel (recall that within our approximations 
the S'-wave pseudoscalar (r/c) and vector {J/tp) channels 
are degenerate) . A common way to display the euclidean 
correlator at a given temperature uses a normalization to 
a so-called reconstructed correlator, which uses a baseline 
spectral function (e.g. the vacuum one) with the Kernel 
K, at the same temperature as in numerator. 



Rc.{t;T) 



JdEaa,iE,T)K.iT,E,T) 
JdEaa,{E,T,,,)ICiT,E,T) 



(36) 



The idea underlying this ratio is to exhibit the tempera- 
ture effects on the in-medium spectral function, (7a{E, T) 
relative to a baseline spectral function, aa{E,T^(,c)j and 
thus to reduce the effects of systematic uncertainties (e.g. 
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discretization effects in IQCD which distort the high- 
energy part of the spectral functions). As pointed out in 
Ref. [3] the spectral function used in the reconstructed 
correlator can have significant impact on the normaliza- 
tion and shape of Ra{T] T). Here, we always use our cal- 
culated vacuum spectral function as baseline, i.e., ricc=0. 
Let us finally elaborate on the diffusion properties of 
a single heavy (anti-) quark which we evaluate in terms 
of our hcavy-Ught quark T-matrix. Following Ref. [53|, 
one may approximate the Boltzmann equation for the 
HQ distribution function in the QGP by a Fokker-Planck 
equation and extract the pertinent thermal relaxation 
rate (inverse relaxation time) as 



Ic = '^Itq^ \\niA{p) 



(37) 



with the friction coefficient 
1 



A{v] 



16i2TT)^0jQ{p)J ^^ 

d^p' (27r)4 



npiujqiq)) 



d^ 



q' 



wq(p') 



1- 



Y,\M\'S^'\q + p- 



UJq{q') 

q -p) 



p-p' 

p2 



(38) 



The invariant amplitude squared, which is summed over 
color, angular-momentum, spin and light-flavor degrees 
of freedom (and averaged over the dc = 6 initial spin- 
color states of a heavy quark), is calculated in terms of 
S- and P-wave on-shell T-matrices as 



Y: |Mp = ^(. - ml + mlf{s -m% + m^ 

x7V/^d,(|f,,o(s)|'+3|f,,i(s) cos(0,„)n (39) 



where^ 

TaA^) = "l(pcm)^^^ Ta,i{E,Pcm,Pcm) ^-(pcm)^^^ (40) 

with center-of-mass (cm) energy and momentum 

£■ = Vs = UJqiPcm) + WQ(pcm) 
Pcm = ^ y TOq + {m^ - s)2 - 2TO|(m2 -t- s) . 

The color degeneracy factors are given by 

do = 1, dg = 3, dg = 6, ds = 8 . (41) 



In Eq. (j38|) . the distribution functions, np, include up 
(m), down (d) and strange (s) quarks in the thermal heat 
bath with inconnng (outgoing) 3-momentum, q{q'). The 



Due to the slightly different definition of the relativistic factors 
in our T-matrix compared to Ref. Q the connection to the cross 
section is modified [53. 



in- and outgoing HQ 3-momenta are p and p' . As an 
extension to previous work [6|, |55| we here distinguish ex- 
plicitly between light- and strange-quark contributions 
(instead of using light quarks with an effective degener- 
acy oi Nf = 2.5). In accordance with HQ spin symmetry 
(as adopted in the quarkonium sector) we assume degen- 
eracy of S*- waves (e.g., pseudoscalar D and vector D* 
mesons) and of P- waves (e.g., scalar Dq and axialvector 
Di mesons). In our numerical calculations below we also 
evaluate the contributions from HQ scattering off gluons. 
In this case, a potential picture is less straightforward. 
Therefore, as in previous work [6|, |55| . we account for 
elastic HQ-gluon scattering via the leading order pertur- 
bative diagrams (including a Debye screening mass) with 
a rather large coupling constant of as = 0.4. 



IV. SPECTRAL FUNCTIONS, CORRELATORS 
AND RELAXATION TIMES 

In this Section we first fix the remaining free param- 
eters, i.e., the bare heavy- as well as light- and strange- 
quark masses and check the resulting spectral functions 
against vacuum spectra of hidden and open heavy-flavor 
mesons (Sec. lIVAp . We then discuss our numerical re- 
sults for quarkonium spectral functions and pertinent eu- 
clidean correlators with emphasis on the uncertainties 
originating from the choice of potential and reduction 
scheme fSec. |IVB|) . Finally, we apply our formalism to 
evaluate HQ thermalization times and diffusion coeffi- 
cients (•Sec. UVCl) . 



A. Vacuum Spectroscopy and Quark Masses 

Let us first focus on the vacuum spectra of charmonia 
and bottomonia to determine the bare masses of charm 
(c) and bottom (6), m° ^, which figure into the expres- 
sion for the effective mass, Eq. (fT2)) . We do this by re- 
quiring the S-wave charmonium (bottomonium) ground 
state to occur at the average mass of r]c and J/tp at 







BbS scheme 


Th scheme 


Potential 1 




1.355 GeV 
4.712 GeV 


1.264 GeV 
4.662 GeV 


Potential 2 




1.402 GeV 
4.768 GeV 


1.293 GeV 
4.700 GeV 



m, = 0.4 GeV, m, = 0.55 GeV 



TABLE I: Dependence of the bare HQ masses on the 3-D 
reduction scheme and underlying lattice potential when ad- 
justing the ground-state quarkonium mass to the average ex- 
perimental value. Also quoted are the effective light- and 
strange-quark masses which have been adjusted to the aver- 
age heavy-light meson ground-state mass. 
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FIG. 8: (Color online) Quarkonium spectroscopy in vacuum for S- and P-wave charmonia (upper panels) and S- and P-wave 
bottomonia (lower panels). 



-3.04 GeV, and at the T(IS') mass of -9.46 GeV (we ne- 
glect hyperfine splittings), see Fig. |S1 Since the entropy 
term in the HQ free energy vanishes in the vacuum there 
is no difference between the free and internal energy. 
The resulting bare-quark masses are compiled in Tab. |l] 
for the two different potentials and reduction schemes. 
They generally fall into the range expected from the bare 
masses quoted by the particle data group jSg] and are also 
consistent with previous T-matrix calculations [7|. The 
spread (in particular the relative one) is somewhat larger 
in the charm sector (ca. 140 MeV) than in the bottom 
sector (ca. 100 MeV), in line with the expectation that 
the 3-D reduction becomes more reliable with increasing 
mass. The mass gap between the ground and first excited 
charmonium state varies rather little between the two 
potentials within a given reduction scheme, (5m^=0.65- 
0.68 GeV (BbS) and 0.54-0.56 GeV (Th). Compared to 
the experimental value of m^i ~ mj/^p ~ 0.59 GeV, the 
Thompson scheme seems to be doing slightly better (the 
bare masses in the BbS scheme tend to be slightly high) . 
The situation is opposite for the pseudoscalar mass split- 
ting between 77c and 77c (25'), where the BbS scheme does 
slightly better (however, the effects of the hyperfine split- 
ting are expected to be larger in the pseudoscalar than 
in the vector case). The Th scheme appears to perform 
somewhat better again for the Xc states, for which the 
BbS scheme overprcdicts the spin-averaged mass by up 



to 0.1 GeV. From these observations we deduce an over- 
all uncertainty of our T-matrix approach of 50-100 MeV 
in the charmonium sector, comparable to corrections one 
expects from hyperfine splittings. 

In the bottomonium sector (lower panels in Fig. [8]), the 
mass gaps between the ground state T and the first ex- 
ited state, as well as between the first and second exited 
state, are reproduced within 30 MeV (BbS) and 70 MeV 
(Th). The differences in the potentials have again only 
minor impact. A similar trend is found for the spin- 
averaged masses of the Xb states: using the BbS scheme 
our results tend to be higher in mass (especially for po- 
tential 2) compared to the experimental values for Jn^^^ 
and 77ij(-t2 , while for the Th scheme we typically obtain re- 
sults 30 MeV below experiment. In addition, for both re- 
duction schems and potentials, we obtain a Xb{3S) state 
right at the continuum threshold. Since there is no exper- 
imental evidence for this state, this could again be indica- 
tive for some over-binding. Recall, however, that we do 
not account for residual B-B interactions in our single- 
channel treatment, which could have a significant impact 
on the spectral function especially close to threshold. As 
to be expected, the difference in the Coulomb term of the 
two vacuum potentials (different as but equal string ten- 
sion) induces larger deviations for the more deeply bound 
bottomonia, while the sensitivity to the reduction scheme 
(static approximation) is reduced. Overall, the accuracy 
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FIG. 9: (Color online) Imaginary part of the charm-light T-Matrix in the vacuum using niq = 0.4 GeV. We show the color 



singlet Qq (upper left plot), 
channels. 



anti-triplet Qq (upper right plot), octet Qq (lower left plot) and sextet Qq (lower right plot) 



of the predictions of our T-matrix approach is at the 10% 
level of the 15-25' mass sphttings. This is of the same 
order (or even below) the observed hyperfine splittings. 
This seems reasonable given that we have neglected both 
spin-spin and spin-orbit interactions at the quark level, 
as well as residual mesonic interactions in DD and BB 
channels. 

Finally, the values of light- and strange-quark mass 
have to be fixed. Since the physics of their effective vac- 
uum masses is rather different than in the HQ sector 
(spontaneous chiral symmetry breaking vs. string break- 
ing), we directly adjust the constituent masses. With 
niq = 0.4 GeV we obtain a 5- wave £)-meson mass of 
2.01(2.02) GeV in the Th (BbS) scheme which coincides 
with the experimental value for the D* meson (but is 
larger than the average D-D* mass by ca. 60-70 McV), 
see Fig. |9l It turns out that both smaller and larger mg 
lead to a larger D-meson mass: in the former case the 
increase in kinetic energy dominates, while in the latter 
case the increase in mass is more important. The re- 
sult for the Z?-meson mass is roughly consistent with the 
string-breaking scale in the HQ potential, in the sense 
that the DD threshold (= twice the ZJ-meson mass) ap- 
proximately coincides with twice the the separation en- 



ergy of the QQ pair plus their bare masses, 

2m.D c^ VirsB) + 2m2 = 2me . (42) 

In this interpretation, the binding energy of the heavy- 
light system should coincide with the constituent light- 
quark mass. This is roughly satisfied in the charm sec- 
tor (m/5 is ^3-10% larger than tUc) while the agreement 
improves in the bottom sector. Interesting effects are 
found in the non-singlet color channels (which will figure 
into our calculations of HQ transport in Sec. IIVCI be- 
low), cf. Fig. ini In the color- antitriplet diquark channel, 
where the Coulomb term brings in half of the attraction 
as in the mesonic (color-singlet) channel, a bound state 
is observed at about niQq ~ 2.1 ± 0.05 GeV, correspond- 
ing to a binding energy of ca. 0.15 GcV. To construct a 
charmed baryon, one may imagine to add another light 
quark with an estimated binding of ~0.25GcV, in anal- 
ogy to the D-uicson. The resulting baryon mass would 
amount to ~2.25GeV, not far from the empirical Ac 
mass of 2.29 GeV. The color-Coulomb is repulsive in the 
sextet and octet channels, implying that the states at 
around ^^2.2 GeV are entirely due to the confining force. 
It is tempting to speculate that the binding of an octet 
and anti-octet (or sextet and anti-sextet), with a bind- 
ing energy comparable to the ground-state charmonium. 
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~0.6GcV, could be a relevant configuration underlying 
the recently discovered X, Y and Z states in the 3.8- 
4.5 GcV mass region. The small widths of these states 
would be naturally explained due to their predominantly 
color non-singlet building blocks, see also Refs. J57H59| . If 
this picture is correct, one predicts further regimes of rich 
spectroscopy for narrow "exotic" 4-quark states around 
masses of 6GeV (2c2c), lOGeV {bbqq, 2b2q, 2b2q) and 
20GeV (2626). 

The empirical heavy-strange mesons, Dg and D*, are 
ca. 100 McV heavier than the non-strange states {D and 
D*). We can reproduce this splitting by choosing a con- 
stituent strange-quark mass of iris = 0.55 GeV, consistent 
with typical values in constituent quark models. Other 
properties of the cs states are quite similar to our find- 
ings for the cq states and will not be reiterated here. This 
also applies to the open-bottom bq and bs states. 



B. Quarkonium Spectral and Correlation Functions 
in the Quark-Gluon-Plasma 

With all parameters fixed and in-medium potentials 
determined we now proceed to compute the spectral func- 
tions of heavy quarkonia in the QGP. These can be tested 
by comparing the pertinent euclidean correlator ratios, 
Eq. (|36|) . to computations of this quantity on the lattice. 
Recent results by Jakovac et al. [601 i^^ quenched QCD 
and by Aarts et al. [61[ for Nj = 2 show small variations 
of about 10% of the correlator ratios for charmonia up to 
temperatures of about 2 Tc, and even less for bottomonia. 
Such a behavior could be semi-quantitatively reproduced 
in several potential model approaches [3, UM, U3l, l63| . 
However, no systematic assessment of relativistic correc- 
tions has been performed in these works. 

We limit our in-medium investigations to the tem- 
perature regime T > 1.2 Tc; closer to T^ the infnitc- 
distance limit of the internal energy, Uoo{T), exhibits 
a rapid increase which is presumably associated with 
the onset of phase-transition dynamics. We do not ex- 
pect pertinent effects to be properly accounted for in our 
current single-channel (QQ) implementation of the T- 
matrix. E.g., close to a second-order phase transition, 
long-range many-body correlations become important, as 
well as new degrees of freedom such as DD channels, 
which are not included here. 

In the following, we divide the presentation into the 
charmonium (Sec lIVBTj) and bottomonium fSec. lIVB2)) 
sectors, followed by a combined evaluation (Sec. lIVB3l) . 



1 . Charmonium 

We begin our in-medium analysis for the S- and P- 
wave channels in the charmonium sector using a small 
"numerical width" of 20 MeV for the c and c quarks (re- 
call the degeneracy of pseudoscalar- vector as well as of 
scalar- axialvector states). Gontrary to the vacuum, we 



now distinguish 2 scenarios depending on whether the 
free [F) or internal (C/) energy is identified as the static 
finite-temperature potential. The results are compiled in 
Figs. HUMII] for U and in Figs. [HH^n] for F as potential. 

Let us first focus on the former case, V{r] T) = 
U{r;T) — Uoo{T). At the level of the in-medium spec- 
tral functions both IQCD inputs and reduction schemes 
share a number of generic trends, all of which were al- 
ready present in the T-matrix calculations of Ref. [7|. 
The S'-wave ground state {rjc, J I '4') survives as a bound 
state up to temperatures of about 2-2.bTc around which 
it merges into the cc continuum. But even at tempera- 
tures as low as 1.2 Tc the medium effects in the potential 
induce a reduction of the binding energy, Eb = 2mc — m^ 
by about a factor of ^2, to Eb — 0.3-0.4 GeV com- 
pared to 0.6-0.8 GeV in the vacuum (for Th and BbS, 
respectively). The effective quark mass at this tempera- 
ture is approximately the same as in vacuum, causing a 
net increase in the mass of the S'-wave ground state to 
TO^(1.2Tc) ~ 3.3-3.4 GeV. For higher temperatures, the 
binding further decreases but this effect is overcompen- 
sated by a reduction in the effective charm-quark mass 
(i.e., in Uoo{T)/2), so that the mass of the state actually 
decreases. Along with the decrease in binding goes a re- 
duction in the strength of the state (= peak height of the 
spectral function at fixed width). The rather subtle dif- 
ferences in the spectral functions become more apparent 
in the euclidean correlator ratios, especially between the 
two reduction schemes (within a given reduction scheme, 
the two different potentials lead to small variations also 
for the correlator ratios). For the BbS scheme, the ra- 
tios deviate by up to 30-40% from one for temperatures of 
1.2-2 Tc- This is too large compared to the 10-15% reduc- 
tion that has been found in IQCD computations |6C)I463 |. 
However, employing the Thompson scheme, the corre- 
lator ratios are within 15% from one, which is better 
in line with IQGD. The technical reason for the differ- 
ence in the correlator ratios between BbS and Thompson 
scheme can be traced back to the larger binding that the 
BbS scheme generates already in the vacuum. This re- 
quires relatively large bare charm-quark masses (recall 
Tab. H]) which in the medium ultimately lead to too large 
a ground-state mass (or continuum threshold) when the 
latter approaches its dissolution (note that in the BbS 
scheme the J/tp (or rjc) mass at 2Tc is still significantly 
above its vacuum mass, while in the Th scheme it has 
dropped below the vacuum value). For the ground-state 
P-wave state (xc) we also find that, right above Tc, it is 
heavier than in vacuum due to the decrease in binding. 
However, due to its relatively the small binding energy 
(in vacuum Eb ^ 0.22-0.25 GeV and 0.3-0.35 GeV in the 
Th- and BbS scheme, respectively) it dissolves just above 
~1.2Tc where it merges into the cc continuum. 

Next we discuss the in-medium charmonium results 
when using F as potential, V{r;T) = F{r;T) - F^{T), 
summarized in Figs. 1121 and [T51 Gompared to the use of 
J7, the in-medium binding is appreciably reduced (recall 
Fig. [S]). For example, the binding energy of the S'-wave 
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FIG. 10: (Color online) Charmonium spectral functions and euclidean correlators in the pseudoscalar (S'-wave) channel at 
various temperatures using the internal energy (U) as potential. Results for two reduction schemes and two different potentials 
are compared. 
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FIG. 11: (Color online) Charmonium spectral functions in the scalar (P-wave) channel at various temperatures using the 
internal energy (f/) as potential. Results for two reduction schemes and two different potentials are compared. 



ground state at 1.2 Tc is reduced by about one order of 
magnitude, and the state dissolves shortly thereafter, at 
~1.3Tc (Fig. [12). The P-wave states have disappeared 
already just above Tc- At the same time the value of the 
potential at infinity provides a smaller selfenergy correc- 
tion (see Fig. [3]) leading to a smaller effective quark mass 
and, consequently, a lowered continuum threshold com- 
pared to using [/. This, in particular, entails no or only 
a small rise in the in- medium mass of the J/?/; above T^. 
For the BbS scheme the drop in effective mass and the re- 
duction in binding nearly compensate each other leading 
to a stable J/i/' mass until dissolution. For the Th scheme 
the smaller bare-quark mass even leads to a net decrease 
of the in-medium J/i/; mass. The euclidean correlator ra- 
tios are again very similar for the different potentials but 
exhibit a significantly different r dependence for the two 
reduction schemes. For the BbS scheme the deviation 
relative to the vacuum correlator is up to ~50% while 
for the Th scheme it is no more than 10%. However, for 
both schemes the temperature evolution is remarkably 
stable, with variations of no more than 15% even in the 
BbS scheme. Thus the rather large overall deviation orig- 
inates from the reconstructed (vacuum) correlator, where 
the problem can be traced back to the large bare-quark 
mass which is required in this scheme due to the large 



vacuum binding energy. 



To further map out uncertainties we consider the in- 
fluence of a quark width on the correlator ratios. In 
Refs. [a, [Ml it has been found that the charm-quark width 
above Tc may be as large as 0.2 GeV. We injected into 
Eqs. mi) a value of Fg = 0.1 GeV for the HQ' width and 
plot the results, using U as potential, in Fig. [TH As an 
immediate consequence, the J/^ width turns out to be 
about twice the single-quark quark width, as to be ex- 
pected. The "dissociation" temperature (loosely defined 
as the temperature where the peak height is reduced to 
less than twice the height of the continuum) decreases sig- 
nificantly compared to the narrow-width approximation, 
to about 1.7 Tc^ the broadening of the resonance struc- 
ture simply accelerates the merging with the continuum 
part. The peak positions (masses) of the narrow-width 
calculation are basically preserved. The correlator ra- 
tios are increased compared to the calculation with small 
quark widths. The magnitude of the effect is relatively 
small for the Th scheme where the spread was already 
small before. For the BbS scheme the increase is more 
significantly: the up to 40% spread in the narrow-width 
calculation is reduced to within 30%. Similar systematics 
are also found when using F as potential. 
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FIG. 12: (Color online) Same as Fig. 1101 but using the free energy (F) as potential. 
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FIG. 13: (Color online) Same as Fig. 1111 but using the free energy {F) as potential. 



2. Bottomomum 

In analogy to the charmonium calculations we supply 
a small "numerical width" of 20MeV for the h and h 
quarks. The in-medium bottomonium spectral functions 
and correlator ratios are compiled in Figs. [T5I+IT61 for U 
and in Figs. [T71+[TS1 for F as potential. 

For the [/-potential, similar to charmonium, the re- 
duction in binding combined with a large effective quark 
mass (similar as in vacuum) leads to an increase in the 
mass of all bottomonium states right above Tc- Within 
the BbS scheme the mass of the lowest T bound state 
varies by less than 100 MeV over the considered temper- 
ature range of 1.2-2 Tc: the lowering of the bb threshold 
and the loss in binding nearly compensate. But even 
at 2Tc a well-defined T(15') bound state persists. The 
T(2S') survives up to a temperature of about l.lTc and 
shows a much larger variation in mass (about 0.5 GeV) 
while the T(3S') basically dissolves at Tc- In the Th 
scheme we observe a similar pattern. For the euclidean 
correlator ratios the calculations within the BbS scheme 
deviate from one by 20-25%, more than seen on the lat- 
tice. However, the relative temperature variations are 
smaller, ca. 10-15%. In the Th scheme the temperature 
variations are further reduced to less than 10%, and also 
the deviations from one are smaller, which is better in line 



with the essentially constant IQCD correlator ratios close 
to one. Further stabilization of our results is conceivable 
with more realistic in-medium widths and/or improve- 
ments in the connection between vacuum and in-medium 
potentials. For the P-wave Xb states the ground state 
melts at about 1.7 Tc while the first excited state dis- 
solves at about 1.2 Tc. 

When using F as potential the reduction in binding is 
again more pronounced, with a dissolution of all excited 
5- wave T's and all Xb states right at Tc- Only the T 
ground state survives until somewhat above 2Tc. Com- 
pared to the calculation with U as potential the strength 
of the state at 2 Tc is reduced by a factor of ~3, indicat- 
ing the lower binding, while its mass is about 200 MeV 
smaller (the loss in HQ mass overcompensates the loss 
in binding energy). Also here the temperature depen- 
dence of the ground-state mass is rather stable. As be- 
fore, the euclidean-correlator ratios are rather sensitive 
to the interplay of HQ mass, quarkonium binding and 
the "polestrength" of the states. The BbS scheme again 
shows appreciable deviations from one for both poten- 
tials, up to 40%, while for the Th scheme these are 10- 
15%. However, the spread in the temperature depen- 
dence is less than 10% for both reduction schemes and 
lattice inputs. We have verified that the inclusion of 
larger HQ widths has effects similar as in the charmo- 
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FIG. 14: (Color online) Same as Fig. 1101 but employing a single-quark width of 100 MeV. 
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FIG. 15: (Color online) Bottomonium spectral functions and euclidean-correlator ratios, using U as potential, for the pseudo- 
scalar channel at various temperatures. We compare BbS- and Th scheme as well as the two different potentials. 



nium case, increasing the correlator ratios by up to 0.1 units at large r. 
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FIG. 16: (Color online) Bottomonium spectral functions using U as potential for the scalar channel at various temperatures 
using the BbS (upper panels) and Th reduction scheme (lower panels). 



3. Discussion of Quarkonium Results 



Let us try to summarize and evaluate the findings in 
the quarkonium sector. Within the Th scheme, all S- 
wave correlator ratios (for both IQCD inputs, for U and 
F, as well as for charmonium and bottomonium) are 
within ca. 15% of one, for all temperatures between 1.2 
and 2Tc. For a given calculation (scenario) the relative 
deviations within this temperature range are, in most 
cases, even smaller, suggesting that the reconstructed 
correlators play a non-negligible role in the absolute un- 
certainty (e.g., "residual" hadronic interactions between 
D and D states in the continuum are not accounted for 
in a single- channel T-matrix as employed here). Within 
the BbS scheme, we generally find larger deviations of 
the correlator ratios from one (by up to ^-^50%); within 
a given scenario, the temperature variations are signif- 
icantly smaller, up to 30% (or even less especially for 
the free energy) . While this may overestimate the uncer- 
tainty associated with the 4D— ^3D reduction scheme (re- 
call that the BbS scheme has a tendency for over-binding, 
even in vacuum, see also the discussion in Appendix |^ , 
it stipulates that the static approximation (especially for 
charmonia) requires further scrutiny if one aims at an ab- 
solute accuracy at the 10% level (applications based on 



the (nonrelativistic) Schrodinger equation are expected 
to be beset with larger uncertainty). We also corrobo- 
rated indications found in Ref. Q that effects of a finite 
spectral width are not negligible either, increasing cor- 
relator ratios at the 5-10% level. Our schematic imple- 
mentation of the in-medium widths has only scratched 
the surface of a full many-body calculation utilizing mi- 
croscopic single-quark spectral functions in the T-matrix 
equation (see. e.g., Ref. |6J| for a recent perturbative 



calculation of the HQ spectral function in the QGP). 

Our analysis corroborates indications from earlier 
studies 0, [ifl [13, [63 that there is currently no deci- 
sive discrimination power between the different scenar- 
ios realized by the use of U ("strong binding") and F 
("weak binding"). When employing U the mechanism 
underlying a constant (or temperature-stable) correlator 
ratio is rather involved, being a combination of 4 com- 
ponents: On the one hand, the binding energies close to 
(but above) Tc are rather large (several 100 MeV), to- 
gether with a large polestrength (due to the steepness of 
the [/-potential at intermediate distances) . On the other 
hand, the effective HQ mass, governed by Uoo{T), and 
thus the QQ threshold energy, are also large (basically 
as in vacuum). With increasing temperature, the binding 
and the polestrength drop, as do the HQ mass and contin- 
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FIG. 17: (Color online) Same as Fig. 115! but using F as potential. 
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uuiTL threshold, thus balancing the (low-energy) strength in the spectral function. On the contrary, when employ- 
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FIG. 18: (Color online) Same as Fig. 1161 but using F as potential. 



ing F, the binding already vanishes just above Tc, and 
the balance in the spectral function upon increasing T is 
between a further loss of strength in the threshold state 
(cusp) and a reduction in the HQ threshold. In par- 
ticular, with the _F-potcntial one does not encounter a 
regime above Tj, with a large variation in binding en- 
ergy. However, going further down in temperature, such 
as regime must inevitably occur when approaching the 
vacuum limit, and similar "complications" as in the U- 
potential calculations above Tc are to be expected. Thus, 
a sensitive test of whether the F-potcntial can be consis- 
tent with IQCD correlators is in the temperature regime 
where the largest variation in binding occurs (which is 
apparently at or below T^). 



Heavy-Quark Diffusion in the Quark-Gluon 
Plasma 



Following the analysis of in-medium quarkonia we now 
turn to evaluating heavy-flavor transport in the QGP. In 
the vacuum we have found that the low-lying D-meson 
spectrum is reasonably well reproduced, but also that 



shallow bound states might occur in colored heavy-light 
two-body channels (recall Fig. IH]). The calculation of the 
heavy-light T-matrix in the QGP requires an additional 
input in terms of the in-medium light quark masses (re- 
call that the in-medium HQ selfenergy is determined by 
the infinite-distance of the free/internal energy accord- 
ing to Eq. ([T^ V Due to chiral symmetry restoration, 
the vacuum constituent-quark mass is expected to ap- 
proach zero; however, the light quarks and gluons most 
likely acquire (chirally symmetric) thermal masses. We 
approximate these by adopting the functional form ex- 
pected from pcrturbative QCD f67]. 



mth(T) 



y^dT, mq{T)^^r, 



'qfl 



'th 



(T) 



mufl=mdfl=0, 771^,0= 0.11 GeV. (43) 



When implemented into a quasiparticlc (QP) description 
of the QGP this form allows to recover an energy den- 
sity, eqp, which is roughly 10-20% below the perturba- 
tive value, independent of temperature [63 |63- We fix 
the strong coupling in Eq. (H51) at g = 2.3, resulting in 
eqp/esB — 0.83, consistent with recent IQCD calcula- 
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FIG. 19: (Color online) Imaginary part of the in-medium on-shell T-Matrix for charm-light quark scattering in the color-singiet 
(upper left), anti-triplet (upper right), octet (lower left) and sextet (lower right) channels. In all cases IQCD potential-1 is used 
for U within the Thompson scheme. Note the factor 100 difference in the y-scales of the upper and lower panels. 



tions [70| for T > 1.4 Tc. This value for g is also compat- 
ible with our pcrturbative calculations for scattering off 
thermal partons {ag ~0.4). 

In Fig. [Tni we compile in-medium charm-light T- 
matrices using U for potcntial-1 within the Th scheme, 
with an in-medium single-quark width of 100 MeV (un- 
certainties due to reduction scheme and potential are 
exhibited in the context of the thermal relaxation rates 
below). In the medium the color-sextet and -octet cor- 
relations fade rapidly due to screening of the attractive 
string part of potential (recall Fig. [5]). The meson (color- 
singlet) and diquark (color anti-triplet) channels feature 
broad "Feshbach resonances" (i.e., resonances at thresh- 
old) up to ^1.5 Tc- Compared to previous T-matrix re- 
sults [y], the diquark state is slightly more robust, due 
to the refined (color-blind) treatment of the string term 
(e.g., the ratio of peak heights for color anti-triplet to 
color singlet at 1.2 Tc is about twice as large). For the 
charm-strange correlations similar patterns are found. 

Next we calculate HQ relaxation rates. The original 
suggestion of nonperturbative effects in HQ diffusion has 
been put forward in Ref. [55| using an effective reso- 
nance model where the masses and coupling strengths 
were free parameters; within reasonable ranges of these, 



a factor 2-4 shorter thermalization times compared to 
pQCD were found. Subsequently, heavy-light T-Matrix 
calculations Q were carried out to render the schematic 
estimates more quantitative (and to check for the ex- 
istence of D-meson resonances in the QGP), roughly 
confirming the results of the resonance model if the U- 
potential is employed. Here, we elaborate for the first 
time a quantitative connection to in-medium quarko- 
nium properties. With the potential and all other pa- 
rameters determined our relaxation rates, A{p), are pre- 
dictions of the approach. They are calculated utilizing 
Eq. ([55]) and displayed in Figs. [50] and [H] as a function 
of the HQ momentum for several temperatures above Tc- 
For completeness, we have added to the T-matrix re- 
sults the contribution from HQ scattering off gluons using 
LO pQCD diagrams (including Debye-screening) with a 
coupling constant as=0.4. At the lowest temperature, 
T = 1.2 Tc, we find 7c = 0.14 - 0.2 fm"!, where most 
of the variation is due to the potential choice while the 
reduction schemes agree within 10% for a given potential 
(pQCD scattering off gluons contributes ca. 0.025 fm~^). 
Thus, in the scattering regime the dependence on the re- 
duction scheme is less pronounced than for bound states 
(see also Appendix |^ . The relaxation rate increases to 
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FIG. 20: (Color online) Charm-quark relaxation rate as a function of 3-momentum calculated in the T-Matrix approach using 
U as potential, compared the LO pQCD with Qs=0.4. A perturbative gluon contribution has been added to the heavy-light 
T-matrix rates. 



0.25-0.33 fm^^ at 2Tc (again with most of the spread 
owing to the difference in the potentials; pQCD scat- 
tering off gluons contributes ca. 0.07 fm"^). The magni- 
tude of the low- momentum relaxation rates at T = 1.2 Tc 
(2 Tc) is a factor 4-5 (2.5-3.5) larger than for a LO pQCD 
calculation for scattering off thermal quarks, antiquarks 
and gluons with as=0.4. They are slightly larger than 
the previous T-matrix results of Ref. [y], where 7c=0.12- 
0.19 fm""'^ has been obtained over the temperature range 
T = 1.1-1.8 Tc for parametrizations of yet two other 
(quenched [Ulel and 7V/=2 [Ulel]) IQCD-based inter- 
nal energies. In the previous calculations [6| a constant 
charm-quark mass of nic = 1.5 GeV was used while we 
here include the in-medium sclfencrgy from the infinite- 
distance limit of the internal (or free) energy. When using 
U, m* is larger than 1.5 GeV up to temperatures of ca. 
1.9 Tc (potcntial-1 with Th scheme, see Tab. U and right 
panel of Fig. [3|). The extra interaction strength in our 
present calculation compared to Ref. |6| is mostly due 
to the color-blind treatment of the string term, particu- 
larly in the diquark channel. We emphasize that the in- 
medium HQ masses as used here are mandatory to main- 
tain consistency with the quarkonium correlator ratios 
where they play a critical role in balancing the changes 



in binding energy. Our investigations actually show that 
the internal energy based on the quenched IQCD input 
from Refs. [l^, ISa] leads to euclidean correlator ratios for 
quarkonia which exhibit a large temperature variation 
(decrease with increasing T) incompatible with IQCD 
results, i.e., well beyond the 30% error margin deduced 
in Sec. IIVB II The large temperature variation (screen- 
ing) in the underlying potential leads to a decrease of the 
thermalization rate with temperature. This feature is 
not confirmed in the more quantitative calculations pre- 
sented here. However, the increase with temperature of 
7c for our T-matrix (plus pQCD gluon scattering) calcu- 
lations is significantly slower than for the LO pQCD cal- 
culations with temperature-dependent Debye mass: for 
T = 1.2 — > 2Tc the former increase by a factor of ~^1.7 
(less for the T-matrix contribution alone), compared to 
a factor of ^2.5 for LO pQCD only (light anti-/quarks 
and gluons). Furthermore, in the T-matrix calculations 
A{p) decreases appreciably with increasing 3-momentum 
while the LO pQCD results are almost constant. This is 
simply due to the fact that with increasing 3-momentum 
the charm quark is less likely to excite a low-energy Fesh- 
bach resonance in collisions with thermal quarks or anti- 
quarks. At high 3-momentum, resummation effects in the 



26 



0.2 



0.1 



Potential 1 
BbS-Scheme 



0.2 



— 0.1 
< 



Potential 2 
BbS-Scheme 



2 3 

p [GeV] 



2 3 

p [GeV] 



0.2 



0.1 



Potential 1 
Th-Scheme 



0.2 



— 0.1 
< 







Potential 2 
Th-Scheme 







1 



2 3 

p [GeV] 

T-Matrix: T=1.2T(; 

pQGD: T=1.2T^ 



4 

T=1 .5 T, 



2 3 

p [GeV] 



T=2T, 



T=1.5T, T=2T, 



FIG. 21: (Color online) Same as Fig. [20] but using F as potential. 
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T-matrix cease and the relaxation rates come closer to 
the LO pQCD results (recall the importance of the proper 
relativistic factors for this behavior). The difference at 
high 3-momentum is mostly due to the smaller value of 
the screening mass of the Coulomb term in our IQCD fit 
relative to the pQCD value, mw^ 



^l + Nj&gT. 



As in Ref. [y], the dominant contribution to the HQ re- 
laxation rate originates from the S-waNe meson (color- 
singlet) meson and diquark (color-triplet) channels, while 
the octet and sextet channels are suppressed (even at 
1.2 Tc), as is immediately inferred from the magnitudes 
of the corresponding T-matriccs in Fig. 1191 The P-wave 
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FIG. 23; (Color online) Bottom-quark relaxation rates as a function of S-momentum calculated in the T-Matrix approach using 
U as potential (plus perturbative scattering off thermal gluons), compared to LO pQCD. 



channels contribute about 30% of the S'-waves. 

When using F instead of U as potential the low- 
momentum charm-quark relaxation rate is reduced by 
approximately a factor of ~2, but still larger by a factor 
of -2 than the LO pQCD results, cf. Fig. [211 Conse- 
quently, they come closer to the LO pQCD results at 
high momentum, even though a significant enhancement 
persists even at p=5 GeV (mostly due to the differences 
in screening mass as mentioned above). 

To put our results in context with other approaches 
we display in Fig. [22] (left panel) the temperature de- 
pendence of the relaxation rate at zero momentum for 
different models. Specifically, we compare our results 
for U and F to LO pQCD, to earlier T-Matrix calcu- 
lations [6| and to estimates from grav ity-gauge duality 
(AdS/CFT) [Zllzi (see also Refs. [lllll for LO calcu- 
lations with running coupling). The uncertainty bands 
associated with our T-matrix calculations are largely gov- 
erned by the differences in the underlying IQCD input. 
As discussed above, the results using U overlap with the 
earlier T-Matrix calculations (where also U has been used 
as potential), especially when the latter would be cal- 
culated with a color-blind string term. When using F 
the results are closer to, but still significantly above, LO 
pQCD. The AdS/CFT rates are markedly larger than 



any of the T-matrix rates, except for extrapolations close 
toT,. 

In the right panel of Fig. [22] we compile the tempera- 
ture dependence of the spatial diffusion coefficients. 



D.= 



T 



rriclc 



(44) 



for the above discussed approaches. We plot Dg in units 
of the thermal wave length of the medium, l/(27rT), 
which renders it suggestive for a connection to the widely 
discussed ratio of viscosity to entropy-density. E.g., in 
kinetic theory for a weakly interacting gas one has the 
approximate relation 



^ « Itd. 



(45) 



In the strongly coupled limit of the AdS/CFT corre- 
spondence, one finds the same parametric dependence, 
albeit with a different numerical coefficient (the conjec- 
tured lower bound of ij/ s = l/An corresponds to to diffu- 
sion at the thermal wavelength, Dg — 1/{2ttT)). Besides 
the quantitative comparison of the Dg values their T- 
depcndence is of particular interest. It is constant for 
AdS/CFT (which has no no scale other than tempera- 
ture; note that the HQ mass is effectively divided out in 
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FIG. 24: (Color online) Same as Fig. (231 but using F as potential. 



Eq. (04])) and almost constant for LO pQCD and the T- 
niatrix approach with F as potential, decreasing by less 
than 5% and up to 30%, respectively, for T = 2 ^ 1.2 T^- 
The variation is larger, ca. 50%, if U is used as potential. 
The largest variation of more than 50% is found with 
the quenched IQCD input [l^ [G^ for U in the previous 
T-matrix calculations, but, as we indicated above, this T- 
depcndence is incompatible with the small temperature 
variation in the euclidean quarkonium correlator ratios. 
Nevertheless, our current, better constrained T-matrix 
calculations support a decreasing trend when approach- 
ing the "critical" temperature from above, as typical for 
many substances at or in the vicinity of a second-order 
transition. 



In Figs. [23| and [M] we display the relaxation rates for 
bottom quarks for the U- and F-potential, respectively. 
The general trends (and quantitative enhancements over 
LO pQCD) are very similar to the charm case so that an 
analogous discussion applies which we do not reiterate 
here. 



SUMMARY AND CONCLUSIONS 



We have set up a common framework to evaluate prop- 
erties of open and hidden heavy-flavor states in the QGP. 
A thermodynamic T-matrix formalism for heavy quarko- 
nia and heavy-light quark interactions has been combined 
with input potentials estimated from heavy-quark free 
energies computed in lattice QCD. Compared to earlier 
calculations, we have refined this link by utilizing a field- 
theoretic ansatz for an effective in-medium gluon prop- 
agator. This enabled the fits to be carried out at the 
level of the color-average free energy while disentangling 
color- Coulomb and confining interactions and thus gain 
insights into their medium modifications via the tempera- 
ture dependence of the associated fit parameters (screen- 
ing masses and coupling strengths). The T-matrix cal- 
culations further allowed us to identify appropriate rela- 
tivistic corrections to the static potential, including dif- 
ferences between vector and scalar interactions for the 
color-Coulomb and confining parts, respectively. E.g., a 
color-Breit correction naturally emerges for the Coulomb 
term. The relativistic corrections are crucial to estab- 
lish quantitative consistency for high-energy scattering 
between pcrturbative QCD and the T-matrix in Born 
approximation. This connection is a prerequisite for a 
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simultaneous treatment of bound and scattering states, 
which was one of the main objectives of our work. 

The bare masses of the charm and bottom quark have 
been fixed to the (spin-averaged) mass of the quarkonium 
ground states, ijc-J/ip and T, in vacuum. The resuh- 
ing mass splittings for the excited states agree with the 
experimental values within ca. ±10%, which is smaller 
than the effects due to hyperfine interactions which have 
been neglected in this work. The largest source of un- 
certainty turned out to be the static reduction scheme 
underlying the scattering equation, while the 2 consid- 
ered lattice potentials induced smaller variations. We 
also verified that the vacuum D- and i?-meson states are 
reasonably well recovered when using typical values for 
the constituent light- and strange-quark mass. As a by- 
product, we found that the scalar treatment of the confin- 
ing force leads to shallow bound states in the color-sextet 
and -octet channels in vacuum, which might be relevant 
for a rather rich spectroscopy of narrow four-quark states 
as discussed in the recent literature. 

Our finite-temperature calculations have been carried 
out within two scenarios of adopting an in-medium po- 
tential from the lattice results, either the free (F) or in- 
ternal ([/) energy. First, we calculated spectral func- 
tions and pertinent euclidean-correlator ratios for heavy 
quarkonia. We confirmed the earlier found trend that 
for F charmonia dissolve rather close to Tc (T^iss— l-2Tc) 
while for U the J/ip may survive up to 2-2.5 Tc- However, 
both scenarios can lead to almost constant correlator ra- 
tios, and thus to agreement with lattice QCD results for 
this quantity. The reason is a small in-medium HQ mass 
correction when using F, while it is larger for U. As in the 
vacuum, we found significant variations due to the static 
reduction scheme, reflected by deviations of up to ^40% 
in the correlator ratios at a given temperature. How- 
ever, within a given reduction scheme, potential choice 
and lattice input, the relative temperature variation of 
the correlator ratios is usually much smaller. This sug- 
gests that future studies should scrutinize corrections to 
the static approximation, but also the role of the recon- 
structed (vacuum) correlator figuring into the denomi- 
nator of the ratios, especially close to threshold where 
hadronic (DD) correlations could become important. 

For heavy-flavor transport in the QGP, the use of 
U leads to a factor of ^2 smaller thermalization times 
and (spatial) diffusion constant compared to F. This is 
largely due to "Feshbach"-type resonances in meson and 
diquark channels up to 1.3-1.5 T^, but nonperturbative 
rescattering strength persists in the heavy-light T-matrix 
for temperatures beyond 2Tc. Even when using F as 
potentail, these effects lead to up to a factor of 2 faster 
thermalization compared to perturbativc scattering. The 
uncertainty due to the reduction scheme is smaller for 
heavy-quark transport coefficients than for quarkonium 
correlator ratios. The screening effects in the interaction 
generate a significant increase of the spatial diffusion con- 
stant (in units of the thermal wavelength) with tempera- 
ture (especially for U) , suggestive for a minimum toward 



Tc. 

Our analyses suggest that a thermodynamic T-matrix 
approach can be used to establish quantitative relations 
between quarkonium survival and heavy-quark transport 
in the QGP. In particular, we have assessed uncertain- 
ties associated with commonly applied static (potential) 
and nonrelativistic approximations. While the latter are 
mandatory in the scattering regime, the former turned 
out to be on the few-tens of percent level, which is rela- 
tively large for the lattice correlator ratios, but relatively 
small in the context of current estimates for heavy-quark 
diffusion coefficients. A pressing issue remains the addi- 
tional uncertainty in the definition of a finite-temperature 
potential, especially when based on model-independent 
input from thermal lattice QCD. 

Several directions for future investigations emerge from 
our studies. As already mentioned, retardation effects 
and the influence of virtual anti-particle contributions 
need to be addressed, especially in the bound-state 
regime, e.g., by replacing the T-matrix by a Dyson- 
Schwinger formalism at flnite temperature. Such stud- 
ies could also facilitate the treatment of heavy-quark in- 
teractions with thermal gluons beyond the perturbative 
level. A more microscopic treatment of the heavy-quark 
width figuring into the 2-particle propagator of the scat- 
tering equation is desirable and in principle straightfor- 
ward. Additional finite-width effects arise via inelastic 
interaction channels, which can be implemented via cou- 
pled channels into the T-matrix equation. For example, 
giuon radiation is expected to become important for high- 
energy charm-quark scattering and/or quarkonium dis- 
solution, while DD or even magnetic charge- anticharge 
states could improve the description around Tc and ex- 
tend it to temperatures below Tc. Heavy-quark suscepti- 
bilities, or more generally correlators of charm quarks 
with conserved charges (e.g., baryon or strangeness), 
which are computed with good accuracy in thermal lat- 
tice QCD, can be calculated with our T-matrix. Here, the 
presence of broad resonances does not necessarily imply 
large signals in such quantities. Finally, the in-medium 
quarkonium and heavy-quark transport properties should 
be implemented into a comprehensive phenomenological 
analysis of pertinent observables in heavy-ion collisions, 
e.g., via rate equations and/or Langevin simulations in a 
realistic bulk medium evolution. This will provide quan- 
titative tests of the equilibrium results in current and 
future experiments and thus advance our understanding 
of strongly coupled QCD matter at temperatures around 
and above Tc- Work along some of these lines has been 
initiated. 
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Appendix A: Differences in the reduction scheme 

In this appendix we discuss differences between the 
BbS and Th reduction scheme. For this purpose we con- 
centrate on the case of heavy quarkonium bound states. 
From Eqs. ^ it follows that the BbS and Th 2-particle 
propagators differ as 



GTiE.k) 



4^Q(fc) 
E + 2ujQ{k) 



G-(i?,fc) 



(Al) 



due to a different treatment of the left-hand cut (vir- 
tual antiparticle contributions). Both reduction schemes 
should give very similar results for the T-Matrix if 



fik) 



4a;Q(fc) 
E + 2LjQ{k) 



1 



(A2) 



This condition is rather well satisfied in the scattering 
region, i.e., above the 2-particle threshold, where the in- 
tegral is dominated by the pole (unitarity cut) of the 
propagator, E — 2u!Q{k) « 0, which implies E ss 2 WQ(fc). 
However, in the bound-state regime, i.e., below thresh- 
old, the situation can be different. For example, in the 
extreme case of i? — >■ the difference between the prop- 
agators becomes as large as a factor of 2, entailing large 
discrepancies in the results for the T-Matrix. Let us try 
to asses the differences more quantitatively for the case 
at hand, i.e., for the binding energy of the charmonium 
ground state in vacuum. Our results for the BbS and 
Th scheme show a ca. 25% difference in the J/^ binding 
energy (the explicit values are quoted in the legend of 
Fig. [26l). As a rough guideline, the influence of G on the 
binding may be estimated by formally writing the solu- 
tion of the T-matrix as T = V/{1 - GV). At the bound- 
state energy, one has GV — 1, and thus a 25% change 



in G approximately "mimics" a 25% stronger potential, 
or binding energy (for the same static input potential). 
Thus, for the BbS propagator in the T-matrix integral 
one should expect, on average. 



fik) = 



4u^Q(fc) 
E + 2u}Qik) 



1.25 



(A3) 



To estimate the relevance of the integration momenta we 
apply a cutoff. A, in the T-matrix equation (|15p and study 
the dependence of /(A) and the J/tJj binding energy on 
this cutoff, as displayed in Fig. [25l Taking as an ap- 
proximative representative momentum the one by which 
half of the binding is built up (A ~ 1 GeV) and evaluat- 
ing the "BbS factor" at this value, one finds /(A) ~ 1.2. 
The magnitude of the deviations between BbS and Th for 
bound states can thus be roughly accounted for and is ex- 
pected to become larger with increasing ratio of binding 
energy to the mass of the constituents. 



Appendix B: Binding energies 

In this appendix we compile the temperature depen- 
dence of the binding energies of the J/tp and T ground 
states. We define the binding energy as the difference be- 
tween the quark-antiquark threshold, 2mQ, and the mass 
of the state in question. In the vacuum the J/^p [rjc) 
binding energy is about 0.65(0.85) GeV for the Th (BbS) 
scheme. The in-medium binding energies are shown in 
Fig. [26] when using U as potential (for F the state al- 
ready dissolves at about 1.3 Tc). One observes that the 
BbS scheme leads to a steeper dependence of the binding 
on temperature compared to the Th scheme while the 
melting temperature is quite similar in both cases. 

A similar pattern occurs for the T ground state, dis- 
played in the left and right panel of Fig. [27] when using U 
and F, respectively Note, however, that the scheme de- 
pendence of the binding energy is significantly reduced in 
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FIG. 27: (Color online) Temperature dependence of the binding energy (es) of the T (or rib) using U as potential (left plot) 
and F as potential (right plot), for various combinations of reduction scheme and IQCD input. 
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FIG. 28: (Color online) Comparison of the pseudoscalar spectral functions for charmonium (upper panels) and bottomonium 
(lower panels) using the full potential (left column) and the color-Coulomb term only (right column). In all cases U has been 
used as potential together with the Th reduction. 



the bottomonium case, to about 10%, reflecting a better 
accuracy of the static approximation due to the larger 
bottom-quark mass, as expected. With the weaker in- 
teraction implicit in F the binding is reduced by about 
a factor of 4. The uncertainty induced by the different 
IQCD inputs is significantly larger than the one caused 
by the reduction scheme. 



Appendix C: Influence of the confining force 

In this appendix we assess the relevance of the confin- 
ing interaction for bound-state formation. Recalling the 
definition of the free energy from Eq. (JT)) , 

Fa{r,T) = Kf(r,T) + F^(r,r) + 2EQ(r) , 
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we repeat our calculations with the string term, V^ , 
switched off while all other parameters are kept fixed. 
The corresponding results for the internal energy, U, arc 
obtained using Eq. ((TU]) with the modified free energy 
(we keep, however, the self-energy from the full calcu- 
lation). The results using U as the potential are pre- 
sented in Fig. [28] using the Thompson reduction. For 
J/'ip (j]c) states the most striking difference occurs in 
the vacuum where without the confining interaction no 
excited bound states are supported and only a modest 
threshold enhancement remains for the ground state. In 
the medium the relevance of the string term gradually 
decreases until the results become similar to the full cal- 
culation for a temperature close to 2 Tc (even though the 
peak height is still smaller). This follows from the sig- 
nificantly stronger screening of the confining relative to 
the Coulomb term (m£)(r) is much larger than mulT) 
above Tc, recall Fig. [2|). Close to Tc half of the binding 
of the J/^p is still supplied by remnants of the confin- 
ing force. These systematics suggest that charmonia are 
rather sensitive to medium effects on the confining force 
in the temperature regime of 1-2 Tc. At first glance it 
might surprise that the calculation without string term 
produces more binding in the medium than in the vac- 
uum. The reason is that, without the string term, the 
internal energy, as given by Eq. (jlOp , leads to a more at- 
tractive potential in the medium (up to ^1.5 Tc) than in 
the vacuum (note that we are still using the large effective 
mass which, of course, is generated by the large-distance 
limit of the string term) . 

For the more tightly bound bottomonia (T) the sensi- 
tivity to the string term is still appreciable. In the vac- 
uum the ground state is only bound by ca. 100 McV while 
the excited states are unbound. In the medium a similar 
trend as in the charmonium sector is observed, in that 
the significance of the string term ceases as temperature 
increases. 

Generally, our findings clearly demonstrate the impor- 
tance of the confining interaction in both charmonium 
and bottomonium spectroscopy, both in vacuum and in 
medium for temperatures of up to ca. 2 Tc- The use of po- 
tentials developed in a perturbative expansion therefore 
omits important physics in the description of quarkonium 
melting in medium. 
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